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The effect of quantum and thermal fluctuations on the phase diagram of spin-2 BECs is examined. 
They are found to play an important role in the nematic part of the phase diagram, where a mean- 
field treatment of two-body interactions is unable to lift the accidental degeneracy between nematic 
states. Quantum and thermal fluctuations resolve this degeneracy, selecting the uniaxial nematic 
state, for scattering lengths 04 > 02, and the square biaxial nematic state for 04 < a2. Paradoxically, 
the fluctuation induced order is stronger at higher temperatures, for a range of temperatures below 
Tc- For the experimentally relevant cases of spin-2 ^'^Rb and ^'^Na, we argue that such fluctuations 
could successfully compete against other effects like the quadratic Zeeman field, and stabilize the 
uniaxial phase for experimentally realistic conditions. A continuous transition of the Ising type from 
uniaxial to square biaxial order is predicted on raising the magnetic field. These systems present a 
promising experimental opportunity to realize the 'order by disorder' phenomenon. 



The ground-state properties as well as the dynamics 
of spinor Bosc-Einstein condensates has been the subject 
of many experimental studies over the past decade (see, 
for instance, [l], 0, H, 0, Q). The spin degree of free- 
dom introduces a rich variety of spin ordered superfluids 
BSSl ^^'^ the interaction between spin order and super- 
fluidity has several interesting consequences, for example 
topological defects of the magnetic texture that trap vor- 
ticity [9| . The precise nature of the spin order realized 
in the ground state depends on the spin dependent two 
body interaction between atoms, parameterized by the 
scattering length in different total-spin channels. Usu- 
ally, a specification of these scattering lengths uniquely 
fixes the spin configuration of the condensate. An inter- 
esting exception occurs in the case of spin two atoms. 
Here, depending on the relative values of the scattering 
lengths in the total spin 0, 2 and 4 channels (00,02,04), 
one obtains either a ferromagnetic state, with a net spin 
moment, a cyclic (or "tetrahedratic" ) state, which breaks 
time reversal symmetry but does not have a net moment, 
or a nematic state which preserves time reversal invari- 
ance. In contrast to the first two states, the nematic 
state requires specifying an additional parameter for its 
description. To visualize the parameter geometrically, 
the state may be represented by its "reciprocal spinor," 
a configuration of four points. The nematic states have 
the symmetry of a rectangle, and the additional parame- 
ter 77 is related to the aspect ratio of this rectangle. This 
parameter is not fixed by the two body interactions, and 
a manifold of accidentally degenerate states remains at 
this level 0. 

We show here that this accidental degeneracy is lifted 
by superfluid phonons, whose velocities depend on the 
precise nematic state being realized. In the low temper- 
ature regime, the quantum zero point energies of the su- 
perfluid phonons leads to an energy splitting between dif- 
ferent nematic states. This is reminiscent of the Casimir 



force arising from photon zero point energies [llj; here 
the 'force' leads to a splitting of nematic states. At higher 
temperatures, thermal fluctuations lift the degeneracy in 
the free energy via entropic effects. Remarkably, we are 
able to derive closed form analytic expressions for this 
splitting in both these limits, yielding the phase diagram 
shown in Fig. [1] Fluctuations select the two high symme- 
try states, the uniaxial nematic (with the symmetry of a 
line) and the square biaxial nematic (with the symmetry 
of the square) in the parts of the phase diagram shown. 
The latter is an elusive phase with a non-abelian homo- 
topy group, that was long sought after in liquid crystal 
systems, but appears naturally here. Two experimen- 
tally relevant systems, spin-2 *^Rb and ^'^Na, are also 
shown on this phase diagram. They are expected to re- 
alize the uniaxial nematic phase. Bose condensates of 
spin-2 ®^Rb have been realized and for the parame- 
ters of that experiment we find that the splitting between 
nematic states from quantum fluctuations to be of or- 
der 0.3 pK while the free energy splitting from thermal 
fluctuations can be as large as 6 pK. The latter should 
be readily observable in experiments with careful control 
over the quadratic Zeeman splitting which competes with 
this fluctuation mechanism and prefers a square nematic 
state oriented in a plane orthogonal to the field. 

An unusual feature of the fluctuation selection is that 
the order is actually stronger at higher temperatures (suf- 
ficiently far below the transition temperature), in con- 
trast to one's intuitive expectation of thermal effects fa- 
voring disorder. This leads to the prediction that for suf- 
ficiently weak background fields, a transition is expected 
on warming the system in the superfluid state, as fluc- 
tuation induced ordering overpowers the external field. 
These effects, both quantum and thermal, are expected 
to be stronger in the ^■^Na. While such 'order by disorder' 
mechanisms have been widely discussed in the theory of 
frustrated magnetism [l^ Il3l |. we are not aware of an 
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unambiguous manifestation of this effect in experiments. 
The spin-2 nematic condensates seem hke a promising 
system to observe this intriguing effect. 

Mean Field Analysis: In the dilute limit, spin-two 
bosons will interact with the pair-potential V{rx — = 
S{ri - r2)(go^o + 92^2 + giPi) where gp = Anh'^aF/m 
and Vf projects into the total spin F state. Such a pair 
potential leads to the interaction Hamiltonian which can 
be conveniently expressed as: 



Hi 



(1) 

where "0 is a five component vector operator whose com- 
ponent '0m (r) destroys a boson in spin state Fz — m 
at position r. Here, ipt is the time reversal of -0, 
namely ('0^^2i 00? "V'j: V'l)- These parameters are 

related to the original ones by a 
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7 = Uao - .94) - 7(52 - 94). The total Hamil 



^o = /d'r|;j|V0|2-;,(0t^). 

The mean-field phase diagram for such a model is 
shown in Fig. [T] 8]. The various extrema correspond 
to the ferromagnetic, nematic, and tetrahedratic phases. 
We will focus on the region where the nematic phase is 
stabilized. The general nematic state, up to an overall 
phase is 0^ = (^^, 0, cost;, 0, and requires speci- 
fication of the additional parameter 77. As rj is varied, 
the rectangle moves through each of the three planes 
xy, yz, xz. V = ^ corresponds to a uniaxial along 
an axis, while rj = (§ ~ ^)7r corresponds to a square 
in one of the coordinate planes. At this level, the differ- 
ent nematic states are degenerate but as demonstrated 
below, fluctuations will remove this accidental degener- 
acy in the nematic phase, producing the dashed phase 
boundary at 02 = 04 (/? = 0). Along this phase bound- 
ary however, the Hamiltonian possesses S0(5) symmetry 
and hence the various nematic phases are exactly de- 
generate. The mean-field energy density in the nematic 
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phase is Eq/V — — 2(0+7) ' ^^^^"^ ^^'^ chemical potential 
itself is given by /i = (a -I- 7)71.0, where uq is the atom 
density. We now analyze the harmonic fluctuations about 
the different nematic ground states. 

Fluctuations: We follow the standard Bogoliubov the- 
ory of the weakly interacting Bose gas to derive the 
spectrum of fluctuations. Consider decomposing the 
field operator into a dominant c-number piece, the mean 
field expectation value, and a smaller fluctuation piece 
ip(r) = ^/rLQ(l)ri 4- (5^(r). Substituting this in the Hamil- 
tonian we expand to quadratic order in the fluctuations 
(linear terms vanish on choosing the mean field solution) 
and diagonalize the resulting Bogoliubov Hamiltonian 
TLb- This is aided by defining new canonical bosons 
via the linear combinations: 61 = —-^[S'^i + S'^^i), 

b2 ^ -75('^*i - h ^ :75((5*2 - <5*-2), and 
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FIG. 1; The phase diagram for a spin-two spinor condensate. 
The dashed line corresponds to the additional phase bound- 
ary dividing the nematic phase into two phases, a uniaxial 
nematic phase and a square nematic phase. 



tonian is then Ti. — Ho + Hint', where the free part is ■'^ 
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Further, we define 
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The Bogoliubov Hamiltonian assumes a particularly sim- 
ple form in these variables. Defining velocities Vj via 
mu^ — [(2/3 — "f)no — 2/?no cos(27; -I- %^)], or equivalently 
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We see that these modes are completely decoupled 
and readily diagonalized by the Bogoliubov trans- 
formation. Physically, the bj modes are connected 
to spin rotations about the three coordinate axes 
with the spectrum LUj{k) — k^v'j + K^k"^ /Arn^ (j — 

1, 2, 3). The p boson generates phase fiuctuations of 
the condensate and hence has the spectrum 1^4 (fc) = 
k^ (a -\- j)no/m + h?k^ /Am? which involves the com- 
pressibility 71o(q; -I- 7). Finally the q mode corresponds 
to fluctuations of the r] parameter and has the spectrum 
uj5{k) = ky/ {—j)no/'m + t9k'^ /Am^. Note, in the long 
wavelength limit all these modes have a linear disper- 
sion. However, for a generic nematic superfluid, only 
the flrst four are goldstone modes of the broken spin 
and phase symmetry; the fifth will actually acquire a 
gap once the degeneracy between the different rj conflg- 
urations is removed. The contribution of these modes 
to the free energy is readily calculated. However, only 
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FIG. 2: Evolution of the free energy difference between uni- 
axial and square biaxial nematic states as a function of tem- 
perature with *^Rb parameters. 



the dispersion of the bi modes depends on the nematic 
parameter 77 and will be responsible for lifting the ac- 
cidental degeneracy. The contribution from these modes 
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where V is the 



volume, will be discussed below. 

Order by Disorder: Before evaluating the free energy, 
we make some general observations arising from symme- 
try. For this purpose the nematic region of the phase dia- 
gram is best described in terms of the variables (ag — a^), 

the X-axis in Fig.fH and the ratio a — The ne- 

■— ' 7 00—04 

matic region corresponds to (oq — 04) < and \a\ < 1. 
Inspecting the velocities Eq. ([2]) one sees that the degen- 
eracy remains unresolved along the line u = 0, which is 
also the phase boundary for the full Hamiltonian due to 
the enlarged symmetry. The symmetry breaking grows 
on moving away from this axis, but has opposite effects 
on either side of this line. Formally, the transformation 
a —^ —a, rj ^ 77 + 7r/2 leaves the approximate free energy 
invariant. 

At zero temperature, the free energy reduces to the 
contribution to the ground state energy for different ne- 
matic states from the zero point motion of the harmonic 
modes. Remarkably, this energy splitting AE{ri) may be 
evaluated in closed form up to an overall constant [l7|: 
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where it is convenient to define the series of functions: 
g^{r), a) = E^=i[l+c^cos(2r?-|-^)]™. The ground state 
is the one with the lowest combination of velocities, and 
turns out to be the square biaxial nematic state {ij — 7r/2) 
for (7 > 0, and the uniaxial state for ct < (77 = 0). The 
latter is relevant to the experimental systems spin-2 *^Rb 
and ^'^Na, which both are believed to have a ~ —I [l8l |. 

At finite temperatures, thermal fluctuations are much 
more effective at lifting this degeneracy; again the state 



with the smaller combination of velocities and hence a 
larger population of thermal excitations will be entrop- 
ically favored. While this may be evaluated numeri- 
cally, there is a broad range of temperatures T ^ 
m.ayi(mv^) where it greatly simplifies. The lower limit is 
a rather small scale, related to the magnetic energy and 
about InK for ®^Rb in 0. In this limit, the linearly dis- 
persing modes obey the equipartition law, and the lead- 
ing term in the free energy that lifts the degeneracy is 
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On evaluating this expression, thermal fluctuations are 
found to lead to the same ground states as quantum fluc- 
tuations. A plot of the free energy variation with 77 for 
CT « — 1 is in the top curve of figure [3l 

At higher temperatures the depletion of the conden- 
sate becomes important. This may be readily taken into 
account by using the temperature dependent condensate 
density obtained from mean field theory, which for a har- 
monic trap is n{T) — 77,9 [1 — (T/Tc)^]. This temper- 
ature dependence has a significant effect only near Tc 
and leads to a vanishing of the above free-energy split- 
ting on approaching Tc, assuming a continuous transi- 
tion. This temperature dependence of the free energy 
splitting between uniaxial and square biaxial states is 
shown in Fig. [2l for the parameters of ®^Rb in . 

These results may be caricatured by an effective Lan- 
dau theory, where sixth order terms are generated by 
fluctuations. While a full discussion of all such terms is 
undertaken elsewhere , we note here that writing the 
five component order parameter as a symmetric traceless 
matrix x 13 > allows for a 6th order term ^|^tr[x^X*'^] 
in the free energy. This term, when evaluated in the ne- 
matic subspace yields ^|xo|^ cos 677, where xo is the mag- 
nitude of the order parameter. Clearly, if A > 0, then 
the free energy is minimized by the square biaxial state 
77 = 7r/2, 7r/6, 57r/6, while if A < 0, then the free en- 
ergy is minimized by the uniaxial state 7; = 0, tt/S, 27r/3. 
This sixth order term vanishes rapidly on approaching 
the transition, accounting for the near degeneracy of ne- 
matic ordered states near this point. Finally, it may be 
recalled that biaxial nematics are rather difficult to ob- 
tain in liquid crystal systems. Within Landau theory this 
is explained by the presence of a cubic term in the free 
energy _Btr [x^] when the nematic order parameter is a 
real symmetric traceless matrix. This can be shown to 
favor the uniaxial state regardless of the sign of B. How- 
ever, in the spinor condensate, the fact that the spin 2 
field also carries a charge quantum number implies the 
absence of such a term (it is a complex symmetric matrix 
with phase symmetry). Hence biaxial nematics may be 
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AE 


AF 


Be 




0.3 pK 


6 pK 


27 mG 




3 pK 


9 pK 


19 mG 



TABLE I: Free energy splittings at T = 0, T = lOOnK, and 
critical fields, calculated for n — 4 x lO^'^cc"^ - the experi- 
mental conditions in Ref. f^. Sodium has a greater disparity 
in its scattering lengths but has a stronger quadratic Zeeman 
efi'ect, explaining its lower critical field. 



realized, and would naturally occur in the a > part of 
the nematic phase region. In addition to the topological 
defects with non-abelian homotopy for the general biax- 
ial nematic, the square biaxial nematic realized here also 
has half superfluid vortices bound to a spin defect. 

Experimental Prospects: Magnetic fields make the 
phase diagram more complex because the splitting due 
to fluctuations receives competition from the quadratic 
Zeeman energy Hq^ ~ —cB^ J (Pxtjj^ F^^ [23|, which fa- 
vors the square uniaxial state. For weak fields the ne- 
matic axis drops into the xy plane ry = 7r/3, 27r/3 be- 
cause {Fx ~ Q) " i\f^j 0, ^,0, has lower Hg^ than 
\ Fz = ) . On increasing the field, the state continuously 
approaches the square state rj = tt/2. In Q a 50-50 split 
between atoms with Fz = ±2 was observed consistent 
with a square state; we believe this to be due to the rela- 
tively large background field reported in the experiments 
that overcomes the order by disorder effect. To observe 
the latter effect, a lower magnetic field as estimated be- 
low is needed. A possible way to observe the uniaxial 
to square state transition is by measuring the fraction 
of atoms with Fz — along a quantization axis parallel 
to the magnetic field as it drops from 25% for uniaxial 
states to zero for the square state (See Fig. H]). Mag- 
nitudes of fiuctuation induced free energy splittings are 
shown in Table[T]for ^^Rb and ^^Na both at zero tempera- 
ture and at lOOnK (roughly half the critical temperature 
for the density in 2]). The critical field beyond which 
the quadratic Zeeman effect dominates is shown in the 
last column. Given that magnetic fields as low as 50mG 
are readily accessible in current experiments we be- 
lieve these effects are within experimental reach. The ef- 
fect could also possibly be used to access the condensate 
temperature by measuring the critical field. In a trap, an 
important consideration is that the linear part of the spin 
wave dispersions must be above the quantization scale of 
the trap, hujtiup < m msLx{vf). 

In conclusion, the role of fluctuations in selecting the 
ground state spin structure of a spin-2 nematic conden- 
sate were studied. Both quantum and thermal fluc- 
tuations were shown to lead to the phase diagram in 
Fig. [TJ For the experimentally interesting cases of ^^Na 
and *^Rb, these effects compete against the quadratic 
Zeeman term, but can predominate for sufficiently weak 
fields. An Ising transition marks this point, which we 
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FIG. 3: Evolution of the uniaxial nematic into square biaxial- 
nematic in a field. The free-energy, with a constant sub- 
tracted, is plotted in pK per atom for the conditions of Ref 
0, for B = 0, 20, 27mG (top to bottom). Reciprocal spinors 
are also shown. At the minima, they deform continuously 
from rectangles into squares as the field is increased. The 
Ising-like symmetry breaking at the critical field corresponds 
to making a choice between which way to deform the square. 




Magnetic Field (mG) 

FIG. 4; Proportion of Fz = component as magnetic field 
is increased to the critical field, for rubidium at density 4 x 
10"cc"^ and T = lOOnK. 



believe can be accessed readily in future experiments. If 
so, this would be one of the first experimental demon- 
strations of the 'order-by-disorder' effect widely studied 
in the context of frustrated quantum magnetism. 
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